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ABSTRACT 

For the first time, all available pseudo-Schwarzschild potentials are exhaustively 
used to investigate the possibility of shock formation in hydrodynamic, invicid, 
black hole accretion discs. It is shown that a significant region of parameter space 
spanned by important accretion parameters allows shock formation for flow in 
all potentials used in this work. This leads to the conclusion that the standing 
shocks are essential ingredients in accretion discs around non-rotating black holes 
in general. Using a complete general relativistic framework, equations governing 
multi-transonic black hole accretion and wind are also formulated and solved in 
the Schwarzschild metric. Shock solutions for accretion flow in various pseudo 
potentials are then compared with such general relativistic solutions to identify 
which potential is the best approximation of Schwarzschild space-time as far as 
the question of shock formation in black hole accretion discs is concerned. 
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1. Introduction 

The process by which any gravitating, massive, astrophysical object captures its surrounding 
fluid is called accretion. Depending on the rotational energy content of the infalling material, 
accretion flows onto black holes may be broadly classifled into two different categories, i.e., 
non- rotating (spherical) and rotating accretion (accretion discs). If the instantaneous dy- 
namical velocity and local acoustic velocity of the accreting fluid, moving along a space curve 
parameterized by r, are u{r) and a(r) respectively, then the local Mach number M(r) of the 
fluid can be defined as M(r) = The flow will be locally subsonic or supersonic according 
to M(r) < 1 or > 1, i.e., according to u{r) < a{r) or u{r) > a(r). The flow is transonic if at 
any moment it crosses M = 1. This happens when a subsonic to supersonic or supersonic to 
subsonic transition takes place either continuously or discontinuously. The point (s) where 
such crossing takes place continiously is (are) called sonic point (s), and where such crossing 
takes place discontinuously are called shocks or discontinuities. It is generally argued that, in 
order to satisfy the inner boundary conditions imposed by the event horizon, accretion onto 
black holes exhibit transonic properties in general, which further indicates that formation of 
shock waves are possible in astrophysical fluid flows onto galactic and extra-galactic black 
holes. One also expects that shock formation in black hole accretion might be a general 
phenomena because shock waves in rotating and non-rotating flows are convincingly able to 
provide an important and efficient mechanism for conversion of significant amount of the 
gravitational energy (available from deep potential wells created by these massive compact 
accretors) into radiation by randomizing the directed infall motion of the accreting fluid. 
Hence shocks possibly play an important role in governing the overall dynamical and radia- 
tive processes taking place in accreting plasma. Thus the study of steady, stationary shock 
waves produced in black hole accretion has acquired a very important status in recent years 
and it is expected that shocks may be an important ingredient in an accreting black hole 
system in general. 

While the possibility of the formation of a standing spherical shock around compact objects 
was first conceived long ago (Bisnovatyi-Kogan, Zel'Dovich, & Sunyaev 1971), most of the 
works on shock formation in spherical accretion share more or less the same philosophy that 
one should incorporate shock formation to increase the efficiency of directed radial infall in 
order to explain the high luminosity of AGNs and QSOs and to model their broad band 
spectrum (Jones & Ellison 1991). Considerable work has been done in this direction where 
several authors have investigated the formation and dynamics of standing shock in spherical 
accretion (Meszaros & Ostriker 1983, Protheros & Kazanas 1983, Chang & Osttriker 1985, 
Kazanas & Ellision 1986, Babul, Ostriker & Meszaros 1989, Park 1990, 1990a). Ideas and 
formalisms developed in these works have been applied to study related interesting problems 
like entropic-acoustic or various other instabilities in spherical accretion (Foglizzo & Tagger 
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2000, Blondin & Ellison 2001, Lai & Goldreich 2000, Foglizzo 2001, Kovalenko & Eremin 

1998), production of high energy cosmic rays from AGNs (Protheroe & Szabo 1992), study 
of the hadronic model of AGNs (Blondin & Konigl 1987, Contopoulos & Kazanas 1995), high 
energetic emission from relativistic particles in our galactic centre (Markoff, Melia & Sarce- 
vic 1999), explanation of high lithium abundances in the late-type, low-mass companions of 
the soft X-ray transient, (Guessoum & Kazanas 1999), study of accretion powered spherical 
winds emanating from galactic and extra galactic black hole environments (Das 2001). 
With equal (if not more) importance and rigor, the question of shock formation in accretion 
discs around Schwarzchild black holes has been addressed by several authors. While the 
initial works in this direction can be attributed to Pukue (1983), Hawley, Wilson & Smarr 
(1984), Ferrari et al. (1985), Swada et al. (1986) and Spruit (1987), it was Fukue (1987) and 
Chakrabarti and his collaborators (Chakrabarti 1989 (C89 hereafter), 1996 and references 
therein, Abramowicz & Chakrabarti 1990, Chakrabarti & Molteni 1993) who were the first 
to provide the satisfactory semi-analytical or numerical global shock solution for transonic, 
invicid, Keplerian or sub-Keplerian rotating accretion around a Schwarzchild black hole. 
Consequently, their works were further supported and improved by several other indepen- 
dent works (Yang & Kafatos 1995 (YK hereafter), Caditz & Tsuruta 1998, Toth, Keppens 
& Botchev 1998). Because of the inner boundary conditions imposed by the event horizon, 
shocks form in BH accretion discs only if the flow has more than one real physical X type 
sonic point (multi-transonic flow). For a particular set of initial boundary conditions, some 
of the above mentioned works report multiplicity in shock location, but such a degeneracy 
can ultimately be removed by local stability analysis, allowing one to assert that only one 
stable shock location is possible. Hereafter, whenever we will use the word 'shock', it is to 
be understood that we will, in general, always refer only the stable shock location unless 
otherwise mentioned. 

The above mentioned works deserve attention because the shocked flows studied there are 
expected to explain the spectral properties of BH candidates. However, thus far in the as- 

trophysical literature, the theoretical study of steady, standing shock formation in accretion 
discs around non-rotating BHs has suffered from two general limitations. Firstly, the shock 
solutions were obtained either on a case by case basis, or, even when successful attempts 
were made to provide a more complete analysis, the boundary of the parameter space re- 
sponsible for shock formation was obtained only for global variation of the total specific 
energy S (or accretion rate M) and specific angular momentum A of the fiow, and not for 
variations of the polytropic constant 7 of the fiow; rather accretion was always considered to 
be ultra-relativistic ^ which may not always be a realistic assumption. As 7 is expected to 



^By the term 'ultra-relativistic' and 'purely non-relativistic' we mean a flow with 7 = | and 7 = | 
respectively, according to the terminology used in Frank et. al. 1992. 
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have great influence on the radiative properties of the flow in general, we think that ignoring 
the explicit dependence of shock solutions on 7 limits claims of generality. Secondly, except 
for YK, all available so called global shock solutions have been discussed only in the context 
of one particular type of BH potential, namely, the Paczyhski & Wiita (1980) potential ($1 
hereafter). Along with the $1, recent studies (Das & Sarkar 2001; hereafter DS, and refer- 
ences therein) enhance the importance of also considering three other pseudo-Schwarzschild 
BH potentials, one ($2 hereafter) proposed by Nowak & Wagoner (1991), and two others ($3 
and $4 hereafter) due to Artemova, Bjornsson & Novikov (Artemova, Bjornsson & Novikov 
1996, ABN hereafter), in mimicking the complete general relativistic space-time for accre- 
tion around a Schwarzschild black hole. Hence we beheve that being restricted to only one 
specific pseudo-Schwarzschild BH potential docs not guarantee the claimed 'global' nature of 
so called global shock solutions present in the literature, rather one must study the transonic 
disc structure as well as shock formation in all available BH potentials to firmly assert the 
ubiquity of shock formation in multi-transonic accretion disc around a Schwarzschild BH. 
In this context, it is to be mentioned here that YK deserves special importance because 
the shock solution due to YK appears to be the only work available in the literature which 
provides the complete general relativistic description of shock formation exclusively for a 
non-rotating BH. Nevertheless, this work deals only with isothermal accretion but one un- 
derstands that global isothermality in BH accretion is difficult to achieve for realistic fiows 
and more general kind of BH accretion is expected to be governed by polytropic equation of 
state. Also YK does not provide the global parameter space dependence of shock solutions. 
A few authors claim to provide the full general relativistic shock solutions for Schwarzschild 
BHs as a limiting case of their results obtained in Kerr geometry (Chakrabarti 1996a,b; Lu, 
Yu, Yuan & Young 1997, LY3 hereafter). In doing so, a number of assumptions were made, 
some of which, however, may not appear to be fully convincing. For example, either the disc 
was supposed to be in conical equilibrium (LY3), which should not be the case in reality 
because the realistic accretion flow should be in vertical equilibrium (Chakrabarti 1996 and 
references therein); some results valid for isothermal accretion were directly applied to study 
the polytropic accretion in an ad- hoc manner (Chakrabarti 1996a), or some of the Newto- 
nian approximations were not very convincingly combined with complete general relativistic 
equations (Chakrabarti 1996b) which does not strengthen their claim for a full general rel- 
ativistic treatment of shock formation. Hence it is fair to say that although literature on 
general relativistic hydrodynamic BH accretion is well enriched by a number of important 
works (The following is an incomplete list of relevant papers on the subject: Novikov & 
Thorne 1973, NT hereafter, Bardeen & Petterson 1975, Abramowicz, Jaroszynski & Sikora 
1978, Lu 1985,1986, Karas & Mucha 1993, Bjornsson 1995, Riffert & Herold 1995, Ipser 
1996, Pariev 1996, Peitz & Appl 1997, Bao, Wiita & Hadrava 1998, Gammie & Popham 
1998, Gammie 1999), no well accepted complete general relativistic global shock solution 
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exclusively obtained for a hydrodynamic accretion disc around a Schwarzschild BH has yet 
appeared in the literature. 

Motivated by above mentioned limitations encountered by previous works in this field, the 
major aim of our work presented in this paper is to provide a generalized formalism which 
is expected to handle the formation of steady, standing Rankine- Hugoniot shock (RHS) in 
multi-transonic hydrodynamic BH accretion fiow and to identify which region of parameter 
space (spanned by every important accretion parameter, namely, A and 7), will be respon- 
sible for such shock formation for all available pseudo- Schawarzschild BH potentials. We 
would also like to compare the properties of multi transonic accretion in these BH potentials 
with complete general relativistic BH accretion as long as the issue of shock formation is 
concerned. 

Hereafter, we will define the Schwarzschild radius as 

_ 2GMbh 

(where Mbh is the mass of the black hole, G is universal gravitational constant) so that the 
marginally bound circular orbit and the last stable circular orbit take the values 2r^ and 
Svg respectively for a typical Schwarzschild black hole. Also, total mechanical energy per 
unit mass on (sometimes called 'efficiency' e) may be computed as —0.057 for this case. 
Also, we will use a simplified geometric unit throughout this paper where radial distance 
r is scaled in units of r^, radial dynamical velocity u and polytropic sound speed a of the 
flow is scaled in units of c (the velocity of light in vacuum), mass m is scaled in units of 
Mbh and all other derived quantities would be scaled accordingly. Also, for simplicity,we 
will use G — c — 1. In next section, we will briefly describe a few important features of the 
four different pseudo-Schwarzschild 'effective' BH potentials used in this work. In §3, we will 
show how we formulate and solve the equations governing multi-transonic BH accretion in 
these potentials which may have shocks. In §4, we will study multi-transonic BH accretion 
using the full general relativistic frame work and will try to argue (in an indirect, but self 
consistant manner) which potential is expected to be the closest approximation of actual 
general relativistic solutions for which regions of parameter space spanned by S, X and 7, as 
long as one concentrates only on shocked flows. Finally in §5 we will draw our conclusion by 
highlighting some of the possible important impacts of study of shock formation on related 
flelds. 



2. Properties of four pseudo- Schwarzschild BH potentials 



Rigorous investigation of the complete general relativistic multi-transonic BH accretion disc 
structure is extremely complicated. At the same time it is understood that, as relativistic 
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effects play an important role in the regions close to the accreting black hole (where most of 
the gravitational potential energy is released), purely Newtonian gravitational potential (in 
the form $jv = cannot be a realistic choice to describe transonic black hole accretion 
in general. To compromise between the ease of handling of a Newtonian description of 
gravity and the realistic situations described by complicated general relativistic calculations, 
a series of 'modified' Newtonian potentials have been introduced to describe the general 
relativistic effects that are most important for accretion disk structure around Schwarzschild 
and Kerr black holes (see ABN for further discussion). Introduction of such potentials 
allows one to investigate the complicated physical processes taking place in disc accretion 
in a semi-Newtonian framework by avoiding pure general relativistic calculations so that 
most of the features of spacetimc around a compact object are retained and some crucial 
properties of the analogous relativistic solutions of disc structure could be reproduced with 
high accuracy. Hence, those potentials might be designated as 'pseudo-Kerr' or 'pseudo- 
Schwarzschild' potentials, depending on whether they are used to mimic the space time 
around a rapidly rotating or non rotating/ slowly rotating (Kerr parameter a ~ 0) black 
holes respectively. Below we describe four such pseudo Schwarzschild potentials on which we 
will concentrate in this paper. It is important to note that as long as one is not interested 
in astrophysical processes extremely close (within 1 — 2 r^) to a black hole horizon, one may 
safely use the following BH potentials to study accretion on to a Schwarzschild black hole 
with the advantage that use of these potentials would simplify calculations by allowing one 
to use some basic features of fiat geometry (additivity of energy or de-coupling of various 

2 2 

energy components, i.e., thermal (-^), Kinetic (^) or gravitational ($) etc., see subsequent 
discussions) which is not possible for calculations in a purely Schawarzschild metric (see §4). 
Also, one can study more complex many body problems such as accretion from an ensemble 
of companions or overall efficiency of accretion onto an ensemble of black holes in a galaxy 
or for studying numerical hydrodynamic accretion flows around a black hole etc. as simply 
as can be done in a Newtonian framework, but with far better accuracy. So we believe that a 
comparative study of multi-transonic accretion flow as well as shock formation using all these 
potentials might be quite useful in understanding some important features of various shock 
related astrophysical phenomena, at least until one can have a complete and self-consistent 
theory of complete general relativistic shock formation exclusively for a Schwarzschild BH. 
However, one should be careful in using these potentials because none of these potentials 
discussed here are 'exact' in a sense that they are not directly derivable from the Einstein 
equations. These potentials could only be used to obtain more accurate correction terms 
over and above the pure Newtonian results and any 'radically' new results obtained using 
these potentials should be cross-checked very carefully with the exact general relativistic 
theory. 
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Paczynski and Wiita (1980) proposed a pseudo-schwarzschild potential of the form 



which accurately reproduces the positions of and rb and gives the value of efficiency to 
be —0.0625 which is in closest agreement with the value obtained in full general relativistic 
calculations. Also the Keplarian distribution of angular momentum obtained using this 
potential is exactly same as that obtained in pure Schwarzschild geometry. It is worth 
mentioning here that this potential was first introduced to study a thick accretion disc with 
super Eddington Luminosity. Also, it is interesting to note that although it had been thought 
of in terms of disc accretion, $i is spherically symmetric with a scale shift of Vg. 
To analyze the normal modes of accoustic oscillations within a thin accretion disc around 
a compact object (slowly rotating black hole or weakly magnetized neutron star), Nowak 
and Wagoner (1991) approximated some of the dominant relativistic effects of the accreting 
black hole (slowly rotating or nonrotating) via a modified Newtonian potential of the form 



$2 has correct form of as in the Schwarzschild case but is unable to reproduce the value 
of Tfy. This potential has the correct general relativistic value of the angiilar velocity fls at 
Tg. Also it reproduces the radial epicyclic frequency k (for r > r^) close to its value obtained 
from general relativistic calculations, and among all BH potentials, $2 provides the best 
approximation for Qg and k. However, this potential gives the value of efficiency as —0.064 
which is larger than that produced by $1, hence the disc spectrum computed using $2 would 
be more luminous compared to a disc structure studied using $1. 

Considering the fact that the free-fall acceleration plays a very crucial role in Newtonian 
gravity, ABN proposed two different BH potentials to study disc accretion around a non- 
rotating black hole. The first potential proposed by them produces exactly the same value 
of the free-fall acceleration of a test particle at a given value of r as is obtained for a test 
particle at rest with respect to the Schwarzschild reference frame, and is given by 



The second one gives the value of the free fall acceleration that is equal to the value of the 
covariant component of the three dimensional free-fall acceleration vector of a test particle 
that is at rest in the Schwarzschild reference frame and is given by 




2 



(16) 




(Ic) 




{Id) 
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Efficiencies produced by $3 and $4 are —0.081 and —0.078 respectively. The magnitude of 
efficiency produced by $3 being maximum,calculation of disc structure using $3 will give 
the maximum amount of energy dissipation and the corresponding spectrum would be the 
most luminous one. Hereafter we will refer to all these four potentials by $j in general, 
where {i = 1,2,3,4} would correspond to $1 (eqn. 1(a)), $2 (eqn. 1(b)), $3 (eqn. 1(c)) and 
$4 (eqn. 1(d)) respectively. One should notice that while all other $iS have singularity at 
r — Tg, only $2 has a singularity at r = 0. It can be shown that for r > 2rg, while $2 is 
flatter compared to purely Newtonian potential $jv, all other $jS are steeper to $Ar- 
At any radial distance r measured from the accretor, one can define the effective potential 
^l^^ {r) to be the summation of the gravitational potential and the centrifugal potential for 
matter accreting under the influence of ith pseudo potential. ^l^\r) can be expressed as: 

^f{r) = Ur) + ^ (2a) 

where A(r) is the non-constant distance dependent specific angular momentum of accreting 
material. One then easily shows that A(r) may have an upper limit: 



A7(r)=ri^$:(r) (26) 

where $^(r) represents the derivative of $i(r) with respect to r. For weakly viscous or 
inviscid flow, angular momentum can be taken as a constant parameter (A) and eqn. (2a) 
can be approximated as: 

$f^(r) = $.(r) + ^ (2c) 

For general relativistic treatment of accretion, the effective potential can nothe decoupled in 
to its gravitational and centrifugal components. For a Schwarzschild metric of the following 
form: 

ds^ = g^ydx^dx'^ 

= _ _ 1^ de + (1 - ^) ' + (^^' + sin^dd(t>^) 

the world line of the accreting fluid is timelike, and the four velocity of the fluid satisfles the 
normalization condition: 

u^u^" = -1 

where u^{u^) is the contra(co)-variant four velocity of the fluid. The angular velocity VL of 
the fluid can be computed as 

^^ut' ^ Xgtt ^ A (r - 1) 
94,4, 
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where A = — ^ is the specific angular momentum which is conserved for fluid dynamics 



as well as for particle dynamics for invicid flow. The general relativistic effective potential 
$^^(r) {excluding the rest mass) experienced by the fluid accreting on to a Schwarzschild 
BH can be expressed as: 



One can understand that the effective potentials in general relativity cannot be obtained by 
linearly combining its gravitational and rotational contributions because various energies in 
general relativity are combined together to produce non-linearly coupled new terms. 
In Fig 1, we plot $^'^'^(r) (obtained from eq. (2c)) and $^^(r) as a function of r in logarithmic 
scale. The value of A is taken to be 2 in units of 2GM/c. ^^^f curves for different $j are 
marked exclusively in the figure and the curve marked by represents the variation of 
$Q^(r) with r. One can observe that ^{-^-^{r) is in excellent agreement with $^^(r), only 
for a very small value of r (r— >rg), ^l-^-^ starts deviating from ^^Q^{r) and this deviation 
keeps increasing as matter approaches closer and closer to the event horizon. All other 
$^-^'^(r)s approaches to ^Q^{r) at a radial distance (measured from the BH) considerably 
larger compared to the case for <l>f-^(r). If one defines A^^^ {r) to be the measure of the 
deviation of ^^■^■^{r) with $^^(r) at any point r, 



One observes that (r) is always negative for (r), but for other <l>^ (r), it normally 
remains positive for low values of A but may become negative for a very high value of A. If 
|A^-^'^(r)| be the modules or the absolute value of A^-^-^(r), one can also see that, although 
only for a very small range of radial distance very close to the event horizon, Al-^^{r) is 
maximum, for the whole range of distance scale while $i is the best approximation of general 
relativistic space time, $2 is the worst approximation and $4 and $3 are the second and the 
third best approximation as long as the total effective potential experienced by the accreting 
fluid is concerned. It can be shown that \A^-^-^{r)\ nonlinearly anti-correlates with A. The 
reason behind this is understandable. As A decreases, rotational mass as well as its coupling 
term with gravitational mass decreases for general relativistic accretion material while for 
accretion in any $j, centrifugal force becomes weak and gravity dominates; hence deviation 
from general relativistic case will be more prominent because general relativity is basically 
a manifestation of strong gravity close to the compact objects. 

Prom the figure it is clear that for ^Q^{r) as well as for all ^^■^■^{r)s, a peak appears close 
to the horizon. The height of these peaks may roughly be considered as the measure of 
the strength of the centrifugal barrier encountered by the accreting material for respective 
cases. The deleberate use of the word 'roughly' instead of 'exactly' is due to the fact that 





Af^(r) = $f^(r)-$^^i(r) 
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here we are dealing with fluid accretion, and unlike particle dynamics, the distance at which 
the strength of the centrifugal barrier is maximum, is located further away from the peak 
of the effective potential because here the total pressure contains the contribution due to 
fluid or 'ram' pressure also. Naturally the peak height for $^^(r) as well as for ^1-^-^{r)s 
increases with increase of A and the location of this barrier moves away from the BH with 
higher values of angular momentum. If the specific angular momentum of accreting material 
lies between the marginally bound and marginally stable value, an accretion disc is formed. 
For invicid or weakly viscous flow, the higher will be the value of A, the higher will be the 
strength of the centrifugal barrier and the more will be the amount of radial velocity or 
the thermal energy that the accreting material must have to begin with so that it can be 
made to accrete on to the BH. In this connection it is important to observe from the figure 
that accretion under $i(r) will encounter a centrifugal barrier farthest away from the BH 
compared to other $jS. For accretion under all <I>iS except <I>i,the strength of centrifugal 
barrier at a particular distance will be more compared to its value for full general relativistic 
accretion. 



3. Multi-transonic flow in various BH potentials and shock formation 

Following standerd literature, we consider a thin, rotating, axisymmetric, invicid steady flow 
in hydrostatic equilibrium in transverse direction. The assumption of hydrostatic equilibrium 
is justified for a thin flow because for such flows, the infall time scale is expected to exceed 
the local sound crossing time scale in the direction transverse to the flow. The flow is also 
assumed to posses considerably large radial velocity which makes the flow 'advective'. The 
complete solutions of such a system require the dimensionless equations for conserved specific 
energy £ and angular momentum A of the accreting material, the mass conservation equations 
supplied by the transonic conditions at the sonic points and the Rankinc Hugoniot conditions 
at the shock. The local half-thickness, hi{r) of the disc for any $i(r) can be obtained by 
balancing the gravitational force by pressure gradient and can be expressed as: 



For a non- viscous flow obeying the polytropic equation of state p — Kp^ is a measure of 
the specific entropy of the flow), integration of radial momentum equation: 




leads to the following energy conservation equation in steady state: 



. 1 2 «e ^2 



(3a) 
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and the continuity equation: 



d 

dr 



[uprhi{r)] = 



can be integrated to obtain the barion number conservation equation: 



Mir. 



1 3 

7 



(36) 



Following C89, one can define the entropy accretion rate M. — MinK^^'^j^^^ which un- 
dergoes a discontinuous transition at the shock location rgh where local turbulence generates 
entropy to increase Ai for post-shock fiows. For our purpose, explicit expression for A4 can 
be obtained as: 



M 



3 

' f 2 



(3c) 



In Eqs. (3a-3c), the subscript e indicates the values measured on the equatorial plane of 
the disc; however, we will drop e hereafter if no confusion arises in doing so. One can 
simultaneously solve Eqs. (3a - 3c) for any particular $i(r) and for a particular set of values 
of A,7}. Hereafter we will use the notation [Pi] for a set of values of {£, A, 7} for any 

particular 

For a particular value of [Pi], it is now quite straight- forward to derive the space gradient of 
dynamical fiow velocity {^)- for flow in any particular ith BH potential $j(r): 



du 
dr 



+ $;(r)) - 



a 
7+1 



u 



2a2 
u(7+l) 



(4a) 



where represents the derivative of Since the flow is assumed to be smooth everywhere, 
if the denominator of eqn. 4(a) vanishes at any radial distance r, the numerator must 
also vanish there to maintain the continuity of the flow. One therefore arrives at the so 
called 'sonic point (alternately, the 'critical point') conditions' by simultaneously making 
the numerator and denominator of eqn. 4(a) equal zero. The sonic point conditions can be 
expressed as: 



1 + 7 



$-(r) -F7$.(r) / A2-Fr^$.(r 



3$:(r)+r<(r) 



(46) 



where the subscript s indicates that the quantities are to be measured at the sonic point (s). 
For a flxed [Pi], one can solve the following polynomial of r to obtain the sonic point (s) of 
the flow: 



27 



r 



$-(r) + 7$-(r) / X'^ + r^^'iir) 



3$;(r)+r$;'(r) 



0. 



(4c) 
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Similarly, the value of {^)- at its corresponding sonic point (s) can be obtained by solving 
the following equation: 



47 
7 + 1 



dr J 



— - 2u, 



7 



7 + 1 



r 



3 $ 

— I 7^ 

r (r 



f du 
\ dr 



+a. 



27 



ir] 



+ 



6(7-1) /<(r)\ 6(27-1 



(i+7)H*:wy 7(7+i)H*;wy 7^(7+1) 



+$ 



3A^ 

r-4 







(4(i) 



Where the subscript (s, i) indicates that the corresponding quantities for any ith potential 
is being measured at its corresponding sonic point (s) and ^'1' {r) = 

For all $j's, we find a significant region of parameter space spanned by [Pi] which allows the 
multiplicity of sonic points for accretion as well as for wind where two real physical inner 
and outer (with respect to the BH location) X type sonic points Tin and Tout encompass one 
O type unphysical middle sonic point rmid in between. For a particular $j, if Ai [Pi] denotes 
the universal set representing the entire parameter space covering all values of [Pi], and if 
Bi [Pi] represents one particular subset of Ai [Pi] which contains only the particular values of 
[Pi] for which the above mentioned three sonic points are obtained, then Bi [Pi] can further 
be decomposed into two subsets Ci [Pi] and T>i [Pi] such that: 



Ci [Pi] C Bi [Pi] only for M (r^,) > M (w) , 



and 



Vi [Pi] C Bi [Pi] only for M (r,„) < M {rout) , 

then for [Pi] e d [Pi] , we get multi-transonic accretion and for [Pi] e Dj [Pi] one obtains 
multi-transonic wind. In Fig. 2. we plot 

(£„ A,) e [Pi] e Ci [Pi] c B^ [P^] and {S„ A,) e [P^] e V, [P,] c B, [Pi] 

for all $i(r)s (marked in the figure) when 7 = 4/3. While the specific energy S is plotted 
along the Y axis, the specific angular momentum A is plotted along X axis. For $i(r), 
the shaded region PQR represents the parameter space spanned by S and A for which 
three sonic points will form in accretion (PQR=Ci ['Pi]) while the wedge shaped un-shadcd 
region PSR represents the parameter space for which three sonic points are formed in wind 
(PSR='Di [Pi]). Similar kind of parameter space division is shown for other $i(r)s as well. 
A careful analysis of Fig. 2 reveals the fact that, at least for ultra-relativistic flow, no region 
of parameter space common to all $j(r) is found for which [Pi] e Ci [Pi] or [Pi] e T>i [Pi]. 
However, signiflcant region of parameter space is obtained for which [Pi] e Ci [Pi] or [Pi] e 
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Vi [Vi] for ^2ir) and $3(r) and a very small region of such common zone in the parameter 
space is obtained (only for extremely low values of the energy and angular momentum of 
the acceting matter) for ^2(1^), ^3(1^) and $4(r). As the flow approaches to its purely non- 
relativistic limit, i.e., as we make 7 — > 5/3, tendency for such mutual overlap of parameter 
space for ^2(1^), ^3(1^) and $4(r) increases. Nevertheless, $1 still remains 'untouchable' by 
^2{f) and $3(r); only a particular region of parameter space (fairly low energy accretion 
with intermediate value of angular momentum) is commonly shared by $4(r) and $i(r). 
One also observe that if £p°-^ and Af*"^ are the maximum available energy and angular 
momentum of the flow for any $i(r) for which [Pi] G Ci [Pi] or [Pi] e T>i [Pi], one can write: 



and 



CTnax cmax . cmax . craax 

^3 -> ^ > ^2 



ymax ^ ymax ^ ymax ^ ymax 



^3 

The above trend remains unaltered as 7 — > 5/3 and we observe that both £:?""^ and Af*"^ 
non-linearly anti-correlates with 7. 

If shock forms in accretion (in this work we will not study the shock formation in wind), 
then [Pj]s responsible for shock formation must be somewhere from the region for which 
\Pi\ e Cj \Pi\i though not all \P^ G Ci [Pi\ will allow shock transition. Using Eqs. (3a - 3c), 
we combine the three standard Rankine-Hugoniot conditions (Landau & Lifshitz 1959) for 
vertically integrated pressure and density (see Matsumoto et al. 1984) to derive the following 
relation which is valid only at the shock location: 



I i-e 

log^ 



(1 - 7) y^^j ^(k^+th) - 6(1 + - Rcomp)-' + (1 + e)-^ = 0, (5) 

where S(^ki+th) is the total specific thermal plus mechanical energy of the accreting fiuid: 



'■■{ki+th) 



, Rcomp and /3 are the density compression and entropy enhance- 
ment ratio respectively, defined as Rcomp = {p+/p~) and (3 = (^Xi+/Xi-^ respectively; 

= 1 — r*^^^''') and r = PRcomp, "+" and "_" refer to the post- and pre-shock quantities. 
The shock strength Si (ratio of the pre- to post-shock Mach number of the flow) can be 
calculated as: 

Si^Rcomp{l + 0)- (6) 

Eqs. (5) and (6) cannot be solved analytically because they are non-linearly coupled. How- 
ever, we have been able to simultaneously solve Eqs. (3 - 6) using iterative numerical tech- 
niques. We have developed an efficient numerical code which takes [Pi] and $j as its input 
and can calculate rsh along with any sonic or shock quantity as a function of [Pi]. It is to be 
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noted that like the references cited in §1, we also obtain multiplicity in the shock location. 
Following C89, we perform the local stability analysis and find that only one rgh which forms 

in between Tout and Vmid is stable for all $j. 

If [Vi] G J-'i [Pi] C d [Pi] represents the region of parameter space for which multi-transonic 
supersonic flows is expected to encounter a RHS at Vsh, where they become hotter, shock 
compressed and subsonic and will again become supersonic only after passing through r^n 
before ultimately crossing the event horizon, then one can also define [Pi] G Qi [Pi] which is 
complement of J^i [Pi] related to Cj [Pi] so that for: 



the shock location becomes imaginary in Qi [Pi], hence no stable RHS forms in that region; 
rather the shock keeps oscillating back and forth. We anticipate that Qi [Pi] is also an 
important zone which might be responsible for the Quasi-Periodic Oscillation (QPO) of the 
BH candidates (see §5). 

Figure 3 demonstrates few typical flow topologies of the integral curves of motion for ultra- 
relativistic (7 = 4/3) shocked flows in various $i's (indicated in the figure). While the 
distance from the event horizon of the central BH (scaled in the units of and plotted in 
logarithmic scale) is plotted along the X axis, the local Mach number of the flow is plotted 
along the Y axis. One can easily obtain such a set of figures for any 7 (and [Pi]) which 
allow shock formation. For all figures, ABCD represents the transonic accretion passing 
through the outer sonic point Tout (marked as B) if a shock would not form. However, 
as Ai of the flow is higher at the inner sonic point rj„ compared to Ai at Tout, the flow 
must encounter a shock at C (the vertical line CE marked by an arrowhead represents the 
shock transition), becomes subsonic and jumps on the branch EF, which ultimately hits 
the event horizon supersonically after it passes through the inner sonic point r^n, which is 
marked on EF by the small circle with a dot at the center. An "*" in the flgure indicates 
the location of the middle sonic point Vmid- The corresponding values of rj„, r^id. Tout, the 
shock location Vgh, and shock strength Si arc indicated at the top of each flgure, while the 
corresponding values of the total specific energy S and angular momentum A for which the 
solutions are obtained, arc indicated inside each figure. GBH represents the 'self-wind' of 
the flow, which, in the course of its motion away from the BH to inflnity, becomes supersonic 
after passing through Tout at B. Collectively, ABCEF represents the real physical shocked 
accretion which connects inflnity with the event horizon. The overall scheme for obtaining 
the above mentioned integral curves is as follows: 

First we compute rj„, r^j^ and rgut by solving eq. (4c). Then we obtain the dynamical 
velocity gradient of the flow at sonic points by solving eq. (4d). For a chosen Mj„ (scaled in 
the units of the Eddington rate MEdd) , we then compute the local dynamical flow velocity 
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u{r), the local polytropic sound speed a(r), the local radial Mach number M(r), the local 
fluid density p(r) and any other related dynamical or thermodynamic quantities by solving 
the eq. (4a-4d) from the outer sonic point using fourth order Runge-Kutta method. We 
start integrating from r^^ut in two difi^erent directions. Along BH, wc only solve for u{r), 
a{r) and M(r) because shock does not form in subsonic flows. However, integration along 
BCD involves a different procedure. Along BCD, not only we compute u{r), a(r) and 
M(r), but also, at every integration step (with as small a step size as possible), we keep 
checking whether eq. (5) is being satisfied at that point. To do so, at each and every point, 
we start with a suitable initial guess value of Rcomp and Si and performs millions of iteration 
to check whether for any set of [Rcomp, ^i], eq. (5) is satisfied at that point and whether for 
such [Rcomp, Si], the value of (3 obtained from eq. (5) becomes exactly equal to -^r-^; in 
other words, whether the entropy generated at that point (if any) becomes exactly equal to 
the difference between the entropies at the inner and the outer sonic points respectively. If 
such conditions arc satisfied at some particular point (point C in the flgure), we argue that 
the shock forms at that point and we can calculate any pre- and the post-shock dynamical 
and thermodynamic quantities at the shock location (i-G-, at C). Once a shock is formed, 
the flow jumps from its supersonic branch BCD to its subsonic branch EG. We again start 
calculating u{r), a{r) and M(r) and any other related flow quantities by solving eq. (4a) 
using fourth order Runge-Kutta method (with the help of eq. (3a-3c) and eq. (4d)), but 
this time from the inner sonic point rj„ of the flow. 

In Figure 4, we present the J^i[Vi]s for all four $j's ($i ^(a), $2 — ^ (b), $3 — > (c), 
and $4 (d)). The specific energy S, specific angular momentum A and the polytropic 
index 7 of the fiow are plotted along the Z, Y and X axis respectively. Each surface for a 
particular $i(r) is drawn for a particular value of 7. While the first surfaces (which have 
the maximum surface areas) on the {S — 7) plane represent ultra-relativistic accretion 
(7 — 4/3), successive surfaces are also shown for higher values of 7, taking a regular interval 
of A7 = 0.025. It is observed that as the flow approaches to its purely non-relativistic 
limit, the area of the {£ — 7) surfaces responsible for shock formation starts shrinking. We 
find that the shock location correlates with A. This is obvious because the higher will be 
the flow angular momentum, the greater will be the rotational energy content of the flow 
and the higher will be the strength of the centrifugal barrier (which is responsible to break 
the incoming flow by forming a shock) as well as the further will be the location of such 
barrier from the event horizon. However, rgh anti-correlates with S and 7. which means 
that for same £ and A, shock in purely non-relativistic flow will form closer to the event 
horizon compared to the ultra-relativistic flow. We also observe that the shock strength Si 
non- linearly anti-correlates with the shock location rgh, which indicates that the closer the 
shock forms to the BH, the higher is the strength Si and the entropy enhancement ratio 
p. The ultra-relativistic flows are supposed to produce the strongest shocks. The reason 
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behind this is also easy to understand; the closer the shock forms to the event horizon, the 
higher will be the available gravitational potential energy to be released and the higher will 
be the radial advective velocity required to have a more vigorous shock jump. Compared 
to $2 and $3, $1 and $4 allow wider spans of 7 as well as A for shock formation. If Smax-, 
Xmax and 'jmax rcprcscnts the maximum values of the corresponding parameters for which 
shock formation is possible, we obtain Smax i^s) > ^maxi^A) > ^max {^1) > ^^max ($2), 

Amax > >^max {"^a) > ^max {^z) > ^max (^2) and 7^ax ($4) > Imax ($l) > Imax {^'i) > 

lmax{^2), respectively. Also we observe that as more and more the flow approaches its 
purely non-relativistic limit, shock may form for less and less angular momentum. For some 
$i(r)s, even a very small amount of angular momentum (A < 1) allows shock formation, 
which indicates that for purely non-relativistic accretion, shock formation may take place 
even for quasi-spherical flow. 



General Relativistic multi-transonic accretion 



Following the arguments provided by NT and Chakrabarti 1996a, we derive the expressions 
for the conserved total specific energy £' (which includes the rest mass energy) and the 
entropy accretion rate M. as: 



and 



S' 



M = 5.657ur 



7(7-1) 



r- 1 



1.25 



7 - (1 + a2) V 1 - 

a' (7 - 1) 



2(7-1) 



[r' + (1 - 7)] 



0.25 



(7a) 



(76) 



.7-(l + a2)_ 

One can see from eq. (7a) that the total specific energy in this case, can not be decoupled 
into various linearly additive contributions of separate physical origin (i.e., kinetic, thermal, 
rotational or gravitational) as it could be done for flows in any pseudo-potential. 
Following the procedure outlined in previous section, one can derive the dynamical flow 
velocity gradient for general relativistic accretion flow as: 



du 
dr 



2r 



2r3-A2 
r3+A2(l-7) 



2r-l 
2r(r-l) 



2a2 
7+1 



5-7r 



+ 



A2-3r2 



4r(r-l) ^ 4[r3+A2(l-r)] 



2a2 

U(«2_l)(-y+l) 



+ 



1-U2 



(8a) 



from which the sonic point conditions comes out to be: 



7 + 1 



7 



1 


2r3-A2 


2r-l 


2r 


r3+A2(l-7) 


2r(r-l) 



5-7r 



4r(r-l) 4[r»+A2(l-r) 



A2-3r2 



(86) 
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The sonic point (s) could be computed by solving the following equation: 

2 



E'^ [rf + A^(l-r.)] 



1 



l-*(rs,A) 



7 (7 - 1) 
,7-^(rs,A), 



= 



(8c) 



where 



and 



r] (rs, A) = 



l + ^*(rs,A) 



*(rs,A) 



1 

2rs 



2rg-A^ 
rJ+A^(l-7) 



2rs-l 
2rs{rs-l) 



A2-3r? 



4r(r,-l) ^ 4[rf+A2(l-r,)] 

The dynamical flow velocity gradient at the sonic point (s) can be obtained by solving the 
following equation: 



2 (27 - 3ag) 
(7 + 1) {ul-lf 



.A) 



7 - (1 + a 



(du 



7 + 1 



a'C (rs, A) 



2e (rs. A) 



40^3 



7-(l 



7 + 1 



2r, 



Sr? - A2 



(r, - 1) r3 + A2 (1 - r,) 



24r2 - A^ (16r, - 13) 



lOrf 



8r3 



A2 (8r2-13r, + 5) 



= 



(Sd) 



where 



e(rs,A) 



5 — 7r 



+ 



A= 



4r(r-l) ' 4[r3 + A2(l -r)]_ 

We solve eq. (8c) and find that hke flows in various $i(r)s, here also a significant region of 
parameter space allows the multiplicity of sonic points for accretion as well as for wind where 
one O type unphysical middle sonic point is flanked in between two X type real physical 
sonic points rj„ and Tout respectively. In Fig. 5 we show the regions of parameter space 
for which multi-transonic flow is obtained for both accretion and wind. The dimensionless 
conserved total speciflc energy S {excluding the rest mass energy) is plotted along the Y axis 
whereas the speciflc angular momentum A is plotted along the X axis. In region bounded by 
PQR and marked by A, three sonic points are formed in accretion and in region bounded 
by PRS and marked by W, three sonic points are formed in wind. While the figure is drawn 
for ultra-relativistic flows, the corresponding regions of parameter space can be obtained for 
any 7. If Smax be the maximum value of the energy and if Xmax and Xmin be the maximum 
and minimum values of the angular momentum respectively, for which three sonic points 
are formed in accretion for any particular 7, we observe that [Smax, Xmax, Amm] non-linearly 
ant i- correlates with 7. In other words, as the flow approaches its purely non-relativistic 
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limit, the area of the region involved in formation of multi-transonic accretion decreases to 
a lower value. 

In Fig. 6, we show the integral curves of motion for general relativistic accretion of ultra- 
relativistic polytropic fluid. For a particular set of [£^,A,7] shown in the figure, ABCD 
represents the accretion passing through the outer sonic point Tout (marked in the figure 
by B) location of which can be found by solving eq. (8c). EBI represents the self- wind. 
Flow along EFGH passes through the inner sonic point rj„ (marked in the figure by F) and 
encompasses a middle sonic point rmid location of which is shown in the figure using a "*" . 
Like Fig. 3, here also we obtain the complete solution topology by integrating eq. (8a) (with 
the help of eq. (7a- 7b) and eq. (8c)) using fourth order Runge-Kutta method. 
If E and 11 be the shock compression and the entropy enhancement ratio (at the shock 
location) for this case (E = 11 = ^^)) can show that the following equation will be 
satisfied when shock forms: 



where T(— /+) and u{—/+) are the pre-/post-shock temperature and dynamical velocities 
of the fiow respectively. However, it is our limitation in this paper that we have not been 
able to formulate or solve any equation which can be used to calculate the shock location in 
general relativistic accretion onto Schwarzschild BHs. Nevertheless, if shock forms in such 
flow (which is, of course, expected), it is obvious that the set of {£, A) responsible for shock 
formation must belong to the region PQR (= [Pgr] G Cqr [Pgr], see §3) of Fig. 5 because 
shock will form only in multi-transonic accretion. The above argument is useful to compare 
accretion flows in various $i(r)s with general relativistic accretion (at least as long as the 
question of of shock formation in multi-transonic flow is concerned) in the following way: 
Suppose for ultra-relativistic flows, we take the region of parameter space [Pi] G J-i [Pi] 
for any $i(r) used in this paper (see Fig. 4), and then superpose that region with PQR 
of Fig. 5 and study which $i(r) provides the maximum overlap between [Pi] e Ti [Pi] 
and [Pgr] G Cqr [Pgr] ■ That particular BH potential is then considered to be the most 
efficient pseudo-potential in approximating the general relativistic, multi-transonic, shocked 
BH accretion. However, such an 'efficiency test' is not entirely unambiguous because as we 
are yet to figure out the exact [Pgr] G J-^gr [Pgr], there may be some possibility that for 
any $i(r), though [Pi] G J-'i [Pi] will overlap with [Pgr] G Cgr [Pgr], but instead of falling 
onto [Pgr] G J-'gr [Pgr] , it will rather overlap with [Pgr] G Qgr [Pgr] because the exact 
boundary between [Pgr] G J-'gr [Pgr] and [Pgr] G Qgr [Pgr] could not be explored in our 
work. Nevertheless, we beheve that still our arguments for the 'efficiency test' is of some use, 
at least until one can find out the exact shock formation zone for general relativistic flow. 
In Fig. 7, we superpose the Fig. 5 on [Pi] G J^i [Pi] for all different $(r)s (marked in the 
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figure) used in our work. Unlike other [Pi] G [Vi\s, [P3] G ^3 [7^3] are drawn using long- 
dashed lines to show its overlap with [7^2] G ^2 [P2] ■ The figure is drawn for ultra-relativistic 
fiow but can also be drawn for other 7s as well. We observe that while [Pi] G J-'i [Pi] has 
excellent overlap (except at very high energy) with [Pgr] G Cgr [Pgr], 1^0 other [P^ G Ti [Pi\s 
have any overlap with it. This leads to the conclusion that at least for ultra-relativistic flow, 
not only $1 (r) is a very good approximation, rather it is the only BH potential to approximate 
for the general relativistic multi-transonic shocked flow. However, as the flow approaches 
its purely non- relativistic limit, we observe that the area of the ovarlaping zone for $i(r) 
decreases with higher 7 and [Pi] G [Pi] is pushed back to overlap rather with [T'cii] G 
T>GR [Pgr\\ hence unlike ultra-relativistic accretion, $i(r) may not be considered such a good 
approximation for purely non-relativistic flows. Also we flnd that a region of low energy - 
high angular momentum [7^4] G JF4 [7^4] starts overlapping with [Pgr\ G Cgr[Pgr\- So for 
high 7 flows, along with $i(r), $4(r) may also be considered as a plausible approximation 
for general relativistic accretion. Shocked flows in $2('") and $3(r) newer show any overlap 
with [Pgr\ G Cgr [Pgr\ for any value of 7; hence these potentials may not be considered to 
mimick the general relativistic multi-transonic accretion flows. 



5. Concluding RemEirks 

In this paper, wc provide a generalized formalism which can formulate and solve the equa- 
tions governing the advective, multi-transonic, hydrodynamic BH accretion in all available 
pseudo-Schwarzschild potentials, which may contain steady, standing, Rankine-Hugoniot 
kind of shocks. We have also formulated and solved the equations governing multi-transonic, 
complete general relativistic BH accretion and wind in Schwarzschild metric and compared 
our pseudo-Schwarzschild solutions with the general relativistic one. The main conclusions 
of this paper are the following: 

(a) We observe that a signiflcant region of parameter space (spanned by the conserved to- 
tal speciflc energy £^ the speciflc angular momentum A and the polytropic index 7 of the 
flow) allows shock formation for all potentials, which leads to the strong conclusion that 
stable, standing RHS are inevitable ingredients in multi-transonic accretion disks around 
non-rotating BHs. The same kind of conclusion was drawn by previous works in this field 
(see §1) only for ultra-relativistic accretion in Paczynski & Wiita (1980) potential, whereas 
we make this conclusion more general by incorporating all available BH potentials to study 
BH accretion for all possible values of 7. 

(b) As the shock forms at a particular radial distance, it is clear that self-similar solutions 
should not be invoked while studying real physical BH accretion and related phenomena. 

(c) It is sometimes argued that a non-standing oscillating shock may modulate the disc spec- 
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trum in order to explain the dwarf novae outburst (Mausche, Raymond & Mattei 1995) or 
QPO (Hua, Kazanas & Titarchuk 1997). In this context, the region of parameter space, for 
which three sonic points are formed in accretion but still no steady, standing shock is found 
(sec §3), can be considered as quite an important zone because [Vi] G Qi[Vi] may provide 
the relevant parameters responsible for such physical processes. 

(d) As long as the shock formation in ultra-relativistic black hole accretion is concerned, the 
Paczynski & Wiita (1980) potential $i(r) is the only pseudo potential which can mimic the 
solutions of general relativistic accretion disc around non-rotating BHs in a very efficient 
way. However, in the purely non- relativistic limit (7 — *-5/3), along with $i(r), another BH 
potential $4(r) proposed by ABN, is also observed to mimic the general relativistic solutions; 
at least for low energy - high angular momentum flows. However, it is interesting to note 
one important feature of the Paczynski & Wiita potential $i(r); like spherically symmetric 
accretion (see DS), for accretion disc also, $i(r) is observed to be in excellent agrcmcnt with 
solutions for ultra-relativistic flow in pure Schwarzschild metric, however, it starts loosing 
(albeight very slowly) its efficiency in mimicking full general relativistic solution with higher 
values of 7, i.e., as the flow reaches its purely non-relativistic limits; although the exact 
reason behind this is not quite clear to us. 

Hot, dense and exo-entropic post-shock regions in advective accretion disks are used as a pow- 
erful tool in understanding the spectral properties of BH candidates (Shradcr & Titarchuk 
1998, and references therein) and in theoretically explaining a number of diverse phenom- 
ena, including millisecond variability in the X-ray emission from LMXBs and the generation 
mechanism for high frequency QPOs in general (Titarchuk, Lapidus & Muslimov 1998 and 
references therein), high energy emission from central engines of AGNs (Sivron, Caditz & 
Tsuruta 1996), formation of heavier elements in BH accretion discs via non-explosive nu- 
cleosynthesis (Mukhopadhyay & Chakrabarti 2000), formation and dynamics of accretion 
powered galactic and extragalactic jets, quiescent states of X-ray novae systems and out- 
flow induced low luminosity of our galactic centre (Das 2001; Das & Chakrabarti 1999). A 
number of observational evidences are also present which are in close agreement with the 
theoretical predictions obtained from shocked accretion model (Rutledge et al. 1999; Muno, 
Morgan & Remillard 1999; Webb & Malkan 2000; Rao, Yadav & Paul 2000; Smith, Heindl 
& Swank 2001). Thus we believe that our present work may have far reaching consequences 
because of the following reasons: 

1. Our generalized formalism assures that that our model is not just an artifact of a par- 
ticular type of potential only and inclusion of every BH potential allows a substantially 
extended zone of parameter space allowing for the possibility of shock formation. 

2. Of course there are possibilities that in future someone may come up with a pseudo- 
Schwarzschild potential better than $i(r), which will be the best approximation for 
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complete general relativistic investigation of multi-transonic shocked flow. In such 

case, if one already formulates a generalized model for multi-transonic shocked accre- 
tion disc for any arbitrary $(r), exactly what we have done in this paper, then that 
generalized model will be able to readily accommodate that new $(r) without having 
any significant change in the fundamental structure of the formulation and solution 
scheme of the model and we need not have to worry about providing any new scheme 
exclusively valid only for that new potential, if any. 

3. Even if someone can provide a completely satisfactory model for shock formation in 
full general relativistic (Schwarzschild) BH accretion, still the utility of this work may 
not be completely irrelevant. Rigorous investigation of some of the shock related 
phenomena is extremely difficult (if not completely impossible) to study using full 
general relativistic framework. Hence one is expected to always rely on these pseudo- 
potentials because of the ease of handling them. For example, it was shown that 
(see §4) the total energy of the general relativistic accretion flow can not be decoupled 
into its constituent contributions, whereas for any kind of pseudo-potential (see §2), all 
individual energy components are linear under addition. This provides enough freedom 
and ease to simply add any extra component in the expression for energy to introduce 
any new physics in the system (radiative forces or magnetic fields for example), which 
is certainly not possible while dealing with full general relativistic astrophysical flows 
around non-rotating BHs. 

Thus, for above mentioned reasons, we believe that compared to all previous works based 
solely on ultra- relativistic accretion in $i, our model is better equipped for handling various 

shock related phenomena. 

It is noteworthy that the idea of shock formation in advective BH accretion is contested 
by some authors (Narayan, Kato & Honma 1997, and references therein). However, the 
fact that their claim against shock formation is, perhaps, inappropriate for many reasons, 
has been shown (Molteni, Gerardi & Valenza 2001) from energy considerations. One can 
understand that the problem of not flnding shocks hes in the fact that non-shock ADAF 
models are, perhaps, unable to produce multi-transonic flows because only one inner sonic 
point close to the BH is explored by such works. 

One can observe that flows characterized by [Vi] G J-'i [Pi] in our work may contain low 
intrinsic angular momentum for some cases (especially for purely non-relativistic flow in 
some of the BH potentials) However, such weakly rotating flows are expected to be allowed 
by nature for various real physical situations like detached binary systems fed by accretion 
from OB stellar winds (lUarionov & Sunyaev 1975; Liang & Nolan 1984), semi-detached low- 
mass non-magnetic binaries (Bisikalo et al. 1998) and supermassive BHs fed by accretion from 
slowly rotating central stellar clusters (lUarionov 1988; Ho 1999 and references therein). 
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Even twenty-eight years after the discovery of standard accretion disc theory (Shakura & 
Sunyaev 1973), exact modehng of viscous multi-transonic BH accretion, including proper 
heating and cooling mechanisms is still quite an arduous task, and we have not yet fully 
attempted this. However, our preliminary calculations show that the introduction of viscosity 
via a radius dependent power law distribution for angular momentum pushes the shock 
location closer to the BH; details of this work will be discussed elsewhere. 
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Figure Captions 



Fig. 1: Effective BH potentials for general relativistic (*3?^^(r)) as well as for pseudo-general 
relativistic($^''^'^(r)) accretion discs as a function of the distance (measured from the event 
horizon in units or rg) plotted in logarithmic scale. The specific angular momentum is chosen 
to be 2 in geometric units. See text for details. 

Fig. 2: Parameter space division for multi-transonic, ultra-relativistic accretion and wind in 
four different pseudo-Schwarzschild BH potentials, see text for details. 

Fig. 3: Solution topologies for multi-transonic, ultra-rclativistic (7 = 4/3) shocked flows in 
different BH potentials as indicated in the figure. Sec text for details. 

Fig. 4: Region of parameter space responsible for shock formation (jFj [Vi]), for four different 
BH potentials $1 (a), $2 (b), $3 (c), and $4 (d). See text for details. 

Fig. 5: Parameter space division for ultra-relativistic, multi-transonic, accretion and wind 
in general relativity. 

Fig. 6: Integral curves of motion for ultra-relativistic, multi-transonic, black hole accretion 
and corresponding 'self-wind' in Schwarzschild metric. See text for details. 

Fig. 7: Comparison of parameter space producing shocked multi-transonic accretion in var- 
ious BH potentials, with parameter space representing multi-transonic black hole accretion 
and wind in general relativity. The figure is drawn for the ultra-relativistic flow, see text for 
details. 
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